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Abstract 

The Casimir effect, a key observable realization of vacuum fluctuations, is usually taught in 
graduate courses on quantum field theory. The growing importance of Casimir forces in microelec- 
tromechanical systems motivates this subject as a topic for graduate many-body physics courses. 
To this end, we revisit the Casimir effect using methods common in condensed matter physics. 
We recover previously derived results and explore the implications of the analogies implicit in this 
treatment. 
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I. INTRODUCTION 



The Casimir effect results from the interplay of zero-point fluctuations and boundary con- 
ditions, and leads to the attraction between two parallel conducting plates in a vacuum.-^ It 
was one of the first predicted, observable consequences of vacuum fluctuations. Traditionally 
it has been taught in graduate courses on quantum field theory- Because the Casimir force 
scales inversely proportionally to the fourth power of the plate separation a, it is only measur- 
able when a is quite small (micron regime). Recently, the Casimir phenomenon has assumed 
a new importance in the design of nanoscale devices .-^^"^ Generalizations of the Casimir 
effect to critical systems with external constraints continue to be fruitful.— iI2i!iLMii£iiii It 
therefore is timely to include it in graduate courses on condensed matter theory. 

In this paper we revisit the Casimir effect, recovering previously derived 
result 9i,iiLi&iiL22i2i with an approach used to calculate observable properties in finite- 
temperature solids, and can therefore be naturally included in the many-body curriculum,— 
possibly serving as a simple pedagogical example of these techniques. 

II. THE CASIMIR COEFFICIENT 

The Casimir effect results from the effect of boundary conditions on the zero-point fluc- 
tuation modes of the electromagnetic field. We will consider the simplest case of two parallel 
conducting plates. The energy, AS, is the finite difference between the zero-point energies 
with and without the plates,- 1 ^ 1 ^ 1 ^ 1 ^ and the force is then the spatial derivative of AS. 
The component of the electric field parallel to the conducting plates must vanish. There are 
two sets of modes that satisfy this condition: the transverse electric (TE) and transverse 
magnetic (TM) modes where the electric or magnetic field are respectively parallel to the 
plates. — The electric field for the transverse electric field modes is given by 

E TE (x, z) = E m{z x g ± )e^sin (™z) , (n > 0) (1) 

q±,n>0 

where x and z are the co-ordinates parallel and perpendicular to the plates respectively, n is 
an integer, and E^ n is the Fourier amplitude of the fields. There is no n = TE mode. The 
corresponding magnetic field is calculated using Faraday's equations V x E = —dB/dt, or 
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FIG. 1: Schematic of the Casimir effect indicating that the normal modes of the electromagnetic 
field between two conducting plates occupy a discrete set of wavevectors. In the transverse electric 
(TE) modes the electric field lies parallel to the plates and n > 1. In the TM modes the magnetic 
field lies parallel to the plates and n > 0. The modification of the frequencies of the zero-point 
fluctuations by boundary conditions changes the energy of the system, creating a pressure on the 
plates. 

B — — V x E. The magnetic field for the TM field modes is given by 
iuj 

B™(x, z) = B ?n(z x «u)e^- cos (^-z) , (n > 0). (2) 

q±,n>0 

We note that there is one extra n = TM mode. The corresponding electric fields are 

-> -» 1 OE -> (? -> 
computed from Maxwell's displacement current equation Vx5 = ^r - ^ or E = (V x 

<r at iuj 

B). The Fourier modes of these fluctuations thus involve a discrete set of wavevectors, 

<?n = (q±,qzn), (3) 

where q zn = mr/a and n is an integer, leading to a discrete set of normal mode frequencies 
^q ± n = c \fq± + Qzn ( see Fig. [I]), where c is the speed of light. 
The zero-point energy of the fields inside the plates is given by 

Q± q_L,n>0 
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where the first term is the zero point energy of the n = TM mode, and the second term 
counts the zero point energy of the TM and TE modes with n > 0. We may rewrite these 
two terms as a sum over all n, both positive and negative, as follows 

hr 00 1 

£ c = i E £\Ai + ^- ( 5 ) 

n=— 00 q ± 

d 2 q± 

In the continuum limit we let Y^ff ~ * A f , „ , where A is the area of the plates, to obtain 

n, 



^=4 < 6 > 



n=— 00 



The quantity £c/A determined from Eq. is dimensionally of the form [Sc/A] = hc[L~ 3 ]. 
Because a is the only length scale in the system, it follows that the change in the zero-point 
energy must have the form 

The fact that this Casimir energy is sensitive to arbitrary interplate separation, a, is a direct 
consequence of the gaplessness, and thus the scale-free nature of the photon field. 

The traditional calculation of K, in the Casimir energy Eq. is performed using a regu- 
larization procedure enforced by a zeta function.- 1 ^ 1 ^ In this paper we present an alternative 
derivation in which we calculate the necessary sums by exploiting the structure of the Bose 
function and the residue theorem of complex analysis. This approach is central to the Mat- 
subara formalism^ used to study many-body systems at finite temperature. The calculation 
of K, by contour integration has been discussed,— 1 ^ and we will adapt this treatment as an 
example of the Matsubara method 22 in a graduate course in many-body physics. Therefore 
we take a brief diversion to describe the technique in general before applying it specifically 
to the calculation of the Casimir coefficient. 

III. THE MATSUBARA APPROACH 

We begin by noting that the Bose function 

nB ^ = eKz/klr _ x ( 8 ) 
has poles on the imaginary axis (see Fig. 2) at z — iv n , where v n = n2irkBT/h, because 



Next we take 

z = %v n + 5, (10) 
where 5 is small so that z is slightly off the imaginary axis so that 

n B {iv z + S) = eM/k l T _ - « (11) 

from which we see that k B T is the residue at each of the poles z = iv n of n B {z\ so that 

^)=J2tu M —Y ( 12 ) 

^— ' ft (2 — lV n ) 

n v ' 

If we have a function, F(z), that does not have poles on the imaginary axis, we can use 
the residue theorem and Eq. ffl2|) to write 

/ dzF{z)n B (z) = 2m V ^F(w n ), (13) 

where C is a contour that encircles the imaginary axis in a clockwise sense, as shown in 
Fig. [21 Equation (fT3"]) can be rearranged to read 

£^" )= ^iS FW " BW - (14) 

which is a key result in the Matsubara approach used to evaluate sums that emerge in the 
study of many-body systems at finite temperatures.— 

We will now apply Eq. ( 1T4T) to the specific case of the Casimir coefficient. To do so, 
we must identify the summations that we need to calculate. We begin with the zero-point 
energy per unit area (J6j) 

Eq he f d 2 q± 



T-T^m^ + ^ (i5 » 

We are interested only in the change in the zero point energy as a result of the plates. In 
the limit of infinite plate separation (a — > oo) the discrete interval in q zn , Aq zn = tt/cl, 
becomes infinitesimal, and the sum over n in Eq. (|15p can be replaced by an integral Yl n = 
V — ^- = — f dq z . Therefore the 
presence of the plates is given by 



— t~ = ~ f dq z - Therefore the change in the zero-point energy per unit area due to the 

' TT/a 7T 



AAA 

where 



hca I j^L I( q± ,a), (16) 



i(q±,°) = ^Ev^+^- / ^& + ql (17) 



z = 2mknT x n 



n=3 



n=2 



n=l 



n=-l 

n=-2 
) n=-3 



c 



FIG. 2: Contour integration path C used to sum over the Matsubara frequencies. 

By making this subtraction, we remove the ultraviolet divergences in the zero-point energy. 
By using the Matsubara method, we can reexpress the sum in Eq. (jTTj) as 



(18) 



where 



F(z) = sjc 2 q\ - z 2 . (19) 
We associate the discrete wavevectors, q zn , with a "Matsubara frequency" cq zn = z/ n . Then 

'2irk B T c s 



cq zn = cn- = n 

a \ n 



where the effective Casimir temperature is given by 

he 



so that 



2a 
1 2a 



c 



k B T c he 



(20) 



(21) 



(22) 



We note that T c scales inversely with the plate separation (T c ~ l/2a). 
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(a) 



-cq s 



C 



cq± 



(b) 



-cqi 



C 



cq±_ F(u + id) 



C F{lu - id) 
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FIG. 3: (a) Contour integration path C used to calculate ^ n F(i^ n ) in Eq. (|23|) . where 



F(z) = y c 2 q"j_ — z 2 , showing branch cuts in F(z) at z = ±cgj_. (b) Distortion of the contour into 
contour C that wraps around the branch cuts of F(z). The integrand of the integral along the 
branch cuts is the difference F(co + iS) — F(uj — iS) between the value of F(z) above and below the 
branch cut. 



Following the Matsubara approach, the sum in Eq. ([18]) can now be rewritten as a contour 
integral^ around the poles at z = iv n of the Bose function n B (z,(3c) yielding 

dz 



hp, 



2ni 



F(z)n B (z,p c ) 



The second term in I(q±, a), Eq. ffTTj) . corresponds to the a — > oo, or (3c 
first term, and thus we may write 



I(q ± ,a)= lim (± f ^F(z) [n B (z, Pc) - n B (z, (5')] 



(23) 
oo limit of the 

(24) 



The subtraction of the Pc — > oo limit of the integrand in Eq. regulates the overall 
integral at large z, guaranteeing that the integrand around a contour at infinity vanishes. 
This procedure permits us to evaluate the integral by distorting the contour around the 
branch cuts in F(z) that extend from z = ±cg_i_ to infinity, as shown in Fig. Q. We then 
obtain 



-cq± poo 
+ f 

cqs 



^L[F{u + id) - F(lu - id)} [n B (co,p c ) - {P c - oo}] . (25) 
Lin 
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To evaluate the branch cut, we note for F(uj ± id) = a/ (cq±) 2 — (u ± iS) 2 , 



F(u + iS) - F{u - iS) = -2i^cu 2 - (cq±) 2 sgn(cu), 



(26) 



is an odd function of u, which permits us to replace n B {uj) by its odd part n B {oj) + - to 
obtain 



+ I ) —lF(u + i6)-F(u- i 6)] 



oo J cq± 



2tu 



{n B (oj, (3 C ) + -}- {Pc ^oo} 



(27) 



Because the integrand is an even function of a;, we can replace this integral by twice the 
integral over positive u to obtain 

1 



I(q±,a) 





roc 


duj 


c 2 . 


'cq ± 


7i 






roc 


duj 


c 2 . 


fcq ± 


71 



oj 2 — c 2 q\ 



(28) 



The change in zero point energy is then given by 
A£<7 „^9^ /" cPq±duj 



-2h 2 (3 c 



(29) 



where we have made the substitution 2a/ c = %fic- By carrying out the integral over q±, we 
obtain 

AS C h 2 p c f 3 

— = y tib(u,Pc). (30) 

Rescaling the integral in Eq. fl30l by changing variables to 2 = hu)/k B T and replacing 
/?c = 2a/ he, we obtain 



2! 
is 



AS, 



c 



dx- 



x 



7T 2 flC 



(31) 



A 67c 2 h 2 (3 c c 2 J e x -\ 720 a 3 

in numerical agreement with previous deri vat ions. XM^^k The associated force per unit 



area is then 



F dAS c 7T 2 he , „ t - 1 Ar/ 2 

- 1.3 X 10~ 3 — rrN/m 2 



A da 240 a 4 (a/|im) r ''"" ^ 

indicating that measurements of the Casimir force must be performed at plate separations at 
or below the micrometer length scale.-^^ 1 ^ When the two conducting plates are parallel, 
the force is attractive, but it can be repulsive in other situations.—* 2 ^ 
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IV. IMPLICATIONS 



We end with a discussion of the broader implications of our approach. Here we use 
the Bose function n B {oj) to impose a spatial boundary condition on quantum fluctuations. 
Key to this treatment is the interpretation of the discrete q zn vectors as a set of Matsubara 
frequencies cq zn = v n so that we are effectively mapping the spatial frequencies of the Casimir 
effect to a discrete set of temporal frequencies in statistical mechanics. This procedure 
leads to an effective Casimir temperature Eq. (I2TI) . ksTc = hc/2a, where Tc is linked to a 
spatial boundary condition. If we now return to statistical mechanics, our treatment of the 
Casimir effect allows us to illustrate a deep relation between finite temperature and boundary 
conditions in time, not in space. This relation arises from the fact that the Boltzmann factor 
e~P B corresponds to the unitary time evolution operator of quantum mechanics 



U[-ih/3] 



(33) 



evaluated at the imaginary Planck time t = —ih/3 = —ih/kBT. Thus statistical mechanics 
can be formulated as quantum mechanics in imaginary time, where the temporal evolution 
occurs along the imaginary time axis with t = —ir, where r e [0, h/3]. 

To give the reader a flavor for this link between finite temperature and temporal boundary 
conditions, we consider a simple harmonic oscillator described classically by the action 



S r 



t-2 



dt 



m(f)' z mujQ<f) 2 



(34) 



where <f>(t), m and ojq are its amplitude, mass, and angular frequency respectively. The 
passage from classical to quantum mechanics is achieved using Feynman's observation^ 
that the amplitude for the oscillator to follow a path <p(t) is exp[iS c \/h]. If \<pi) is the 
eigenstate of displacement, then the transition amplitude between |0i) and |02) is the sum 



W(*)} 



over all paths that link <f>\ and 02- We pass from real to imaginary time, and identify t 



(35) 



-IT 



and the Euclidean action Se by — S c \ 

n 



— Se- By making the necessary replacements, we 
n 



find that the imaginary time evolution from \<j>-y) to |0 2 ) is associated with the amplitude 

(36) 



^\e- m \<h) 



Yl ex P 

#o)=<fc, <t>(hp)=<fo} 



-^Se 
n 



To obtain the associated partition function, we must take the trace over this matrix 

1 



Z = Tr[e~ pH ] 



d<f>(<f>\e- pji \<l>) = / d(f) 



ex P 

<)>=4>{o)=<l>{np) 



n 



-Sf 



e*p (-j:S B [{<l>}} 



(37a) 
(37b) 



where the sum is over all periodic paths that satisfy 0(0) = 4>(hf3). We see that a finite 
temperature is formally equivalent to a periodic boundary condition in imaginary time. We 
next apply these ideas to our illustrative case. 

In the passage from real to imaginary time, the quantum mechanical amplitude has been 
replaced by a "probability" £>[{</>}] oc e~ SE ^^ h . To calculate this function for our simple 
example, we decompose the displacement in terms of its normal modes (Fig. H|) by Fourier 
transforming in imaginary time 



w 2 ? 



(38) 



where v n = 2iin/fij3 and 0* = 0_ n because 0(r) = 4>{j)* is real. We can then write the 
Euclidean action as a sum of contributions from each normal mode 



Se[{4>}} _ m_ 

h ~ 2h 



(39) 



so that the probability amplitude p[{0}] = Z l e Se ^ h factorizes into a product of Gaussian 
distribution functions for each normal mode: 

|2 



p[{0}] oc exp 



E 



where 



m 2 ) 



n 



(40) 



(41) 



is the temperature-independent variance of each normal mode. 

Let us now see how boundary conditions in time cause the system to become "hot" for 

our example. More specifically consider the variance in the displacement 

-hp 

+,?• (42) 



i r n @ i i 



which has been rewritten in terms of the normal modes so that the average potential energy 
is given by 



(V) 



mujn 



2 > 



1 



(43) 
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FIG. 4: Periodic paths in imaginary time for the harmonic oscillator can be decomposed into their 
Fourier modes, with amplitudes <p n = 4>' n + i<j>'^. 



(a) 



-Uq 



w 



(b) 



-CJ 



FIG. 5: (a) The contour integration path C used to evaluate Eq. (|45p . (b) The distortion of the 
contour to form the path C' around the poles at z = ±a>o- 



We can calculate this expression using contour integration methods. We rewrite it as a 
clockwise contour integral around the poles of n(z) = [e az / fcsT — at z = %v n (Fig. 0(a)) 



(V) 



dz 
c 2ni 
dz 



n B (z 



1 J c ,2 m nB[Z) (uj 2 -z 2 ) 



(44a) 
(44b) 



where C runs clockwise around the poles at z = ±co>o (Fig. 0(b)). The resulting integral is 
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(V) 



n B {uJ ) n B (-uj ) 



2Un 



n B (uj ) + 2 



(45a) 



From the virial theorem^ we expect a similar expression for the kinetic energy so that the 
total energy of the harmonic oscillator is 



(f + V) = hu 



n B M + 2 



(46) 



The second term is due to T = zero-point fluctuations, and the first describes the thermal 
excitations of the oscillator. From this derivation we see that the latter result from the 
redistribution of normal mode energies, and are a direct consequence of the imposition of a 
boundary condition in (imaginary) time. 

Blackbody radiation is a well-known example of thermally excited oscillators. We can 
easily generalize Eq. ()46l) to the electromagnetic vacuum, an ensemble of harmonic oscilla- 
tors, by replacing ujq — > cq, so that the energy density at a finite temperature is given by 



St 
V 



n B (cq) 



1 

2 



d 3 q 
(2^ 



Ticq 



n B (cq) 



(47a) 



where the factor of two is due to the different polarizations. With the substitution x = hcqfi, 
the thermal energy density takes the form 

AS T {k B Tf 



dxx n B (x). 



(48) 



V ~ 7T 2 C 3 ^ 3 

We compare Eq. (l4*8|) with the energy density obtained in the Casimir effect in Eq. (1301) . 
With fie 1 = k B T c = hc/2a, the Casimir energy density becomes 



AS, 



c 



(k B T ( 



c) 



dx x 3 n B (x). 



(49) 



Aa 3n 2 c 3 h 3 

The similarity between the Casimir and the blackbody energy density is testament to their 
common origin as boundary-condition effects. The sign difference is due to the subtle dis- 
tinctions between imaginary and real time, and the consequences for the kinetic energy. As 
an aside, we note that the blackbody radiation pressure, P = A£t/3V, has the same pref- 
actor (apart from the sign) as the Casimir pressure.— Traditionally we think of black body 
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radiation as resulting from an excitation of thermal modes. Our calculation shows that the 
Casimir effect and blackbody radiation are both consequences of boundary conditions and 
the redistribution of zero-point fluctuation modes in the vacuum. 

Recent experiments have observed the Casimir effect between parallel plates with a one 
micron separation (a = 1 jum) .- The corresponding "Casimir temperature" for these 
experiments is 

T C = ^-~ 1000K. (50) 
lokb 

The Casimir effect at these length scales couples to the same photons that predominate in 
the blackbody spectrum at 1000 K. The boundary conditions imposed by the two phenomena 
on the electromagnetic field are almost identical. 

More generally, zero-point fluctuations play a major role at quantum phase 
transitions.— 1 ^ 1 ^ The effect of finite temperature in the vicinity of a (T = 0) quantum 
critical point is the temporal analog of the Casimir phenomenon, a "Casimir effect in time," 
where temperature imposes temporal constraints on critical zero-point fluctuations. As we 
have discussed, there is an intimate connection between a finite temporal dimension and 
a nonzero temperature in a quantum system,— i22i2L22iM and this connection has many ob- 
servable consequences on thermodynamic quantities for quantum critical systems at nonzero 
temperatures.— Heuristically, this link between temperature and a boundary condition in 
time can be understood within the framework of the Heisenberg uncertainty principle 

A* ~ JL, (51) 

where a thermal energy fluctuation leads to an upper cutoff in time, the Planck time, that 
is inversely proportional to the temperature. More formally, as we have seen in our simple 
example, finite-temperature emerges in a path integral framework as a periodic boundary 
effect in imaginary time, which becomes particularly important near a quantum critical 
point where there exist quantum fluctuations on all spatial and temporal scales. Here finite- 
temperature corresponds to the redistribution of quantum zero-point fluctuations due to 
the imposition of external constraints, and thus is analogous to the Casimir effect for two 
parallel metallic plates in vacuum. Running this argument the other way, we note that the 
removal of temporal modes by periodic finite boundary conditions generates a temperature 
(and thus entropy and thermal energy) in a system near a quantum critical point. As an 
aside, we remark that finite-temperature effects resulting from boundary constraints have 
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been discussed in the context of astrophysics where blackbody radiation and event horizons 
have been linked via the Unruh effect.—^ 

In conclusion, we have revisited the Casimir effect with an approach used to study con- 
densed matter systems at finite temperature. We recovered results previously derived by 
other methods and also discussed the physical implications of analogies implicit in this treat- 
ment. Our hope is that this presentation will make this topic straightforward to include in 
a graduate many-body course for future condensed matter physicists. 
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